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The interpretation of extensive air shower (EAS) measurements is strongly dependent on the hadronic interaction models 
used for simulating reference showers. We study the importance of low-energy hadronic interactions in simulated air showers, 
generated with the simulation package CORSIKA, for the observed characteristics of extensive air showers. In particular we 
investigate in detail the energy and the phase space regions of secondary particle production which are most important for muon 
production. This phase space region is covered by fixed target experiments at CERN. In the second part of this work we present 
preliminary momentum spectra of secondary tt + and ir~ in p+C collisions at 12 GeV/c measured with the HARP spectrometer 
at the PS accelerator at CERN. In addition we use the new p+C NA49 data at 158GeV/c to check the reliability of hadronic 
interaction models for muon production in EAS. Finally, possibilities to measure relevant quantities of hadron production in 
existing and planned accelerator experiments are discussed. 



1. Relation of muons in EAS to hadronic interac- 
tions 

In order to extract information (energy and particle 
type) on the primary particle from ground based air 
shower measurements, simulation of air showers us- 
ing electromagnetic and hadronic interaction models 
are necessary. These simulations show uncertainties 
which come mainly from hadronic interaction mod- 
els. Muons measured at ground level are one of the 
main ingredients to infer energy and mass of the pri- 
mary particle of an air shower. Furthermore muons 
are sensitive to the characteristics of hadronic interac- 
tions. 

One aim of this work is to specify the low-energy 
hadronic interactions which are important for the 
muon production in EAS. We simulate extensive air 
showers with a modified version of the simulation 
package CORSIKA (TJ. This special version stores 
the information of the last hadronic interaction where 
mesons are produced which decay into muons. The 
initiator hadron of the last hadronic interaction is 
called the grandmother particle and the meson which 
decays into a muon is named mother particle. For the 
following analysis vertical proton showers with a pri- 
mary energy of 10 15 eV are used. In this CORSIKA 
simulation, the hadronic interaction model GHEISHA 
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[2 1 is applied for energies up to 80 GeV and QGSJET- 
01 [ 3 1 for higher energies. Simulations with other 
models are in preparation. In Fig.[T]the energy spec- 
tra of grandmother particles of different particle types 
are shown. Most of the grandmother particles are 
pions but also about 20% are nucleons and only a 
few are kaons. The most important energy region 
reaches from 8 GeV up to 1000 GeV. The transition 
from the low to the high-energy model is visible in 
all spectra as a step at 80 GeV. The largest fraction of 
mother particles are pions with around 90%, the rest 
are mostly kaons. For more details see j4j. 

2. Comparison: EAS and fixed target experiment 

In particular the forward hemisphere of the phase 
space of secondary particles is important for muon 
production in EAS as shown in detail in chapter [3] 
Therefore fixed target experiments are well suited to 
measure reactions comparable to the hadronic inter- 
actions in EAS. In contrast to collider experiments, in 
fixed target experiments the forward direction of the 
phase space of secondary particles is accessible and 
nuclear targets can be used. Comparing the last inter- 
action in EAS with collisions studied at accelerators, 
one has to keep in mind that the grandmother particle 
corresponds to the beam particle and the mother par- 
ticle is equivalent to a secondary particle produced in 
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Figure 1 . Energy distribution of grandmother parti- 
cles in vertical proton showers with a primary energy 
of 10 15 eV simulated with CORSIKA (GHEISHA 
and QGSJET-01 implemented) 0). 

e.g. a minimum bias p-N interaction. The most prob- 
able energy of the grandmother particle is within the 
range of beam energies of fixed target experiments 
e.g. at the PS and SPS accelerators at CERN. 

3. Relevant phase space regions 

Motivated by the availability of protons as beam 
particles at accelerators we consider only those last 
interactions in EAS that are initiated by nucleons. The 
phase space regions of mother particles produced in 
the last interaction in EAS are shown in Fig. [2] The 
relevant phase space is given as box histogram. The 
size of the boxes indicates the relative importance of 
the phase space region. The momentum of the grand- 
mother particle is given as abscissa and the mother 
particle observable as ordinate. In both figures, ex- 
isting fixed target p+C data are indicated by shaded 
(colored) regions. In the past the only measurement 
of p+C collisions, which was not limited to a fixed 
angle, was the experiment done by Barton et al. 0. 
These data were collected using the Fermilab Single 
Arm Spectrometer facility in the M6E beam line. A 
proton beam with a beam momentum of lOOGeV/c 
and a thin carbon target (1.37 g cm -2 ) was used. The 
phase space of the secondary particles (pions, kaons, 
protons) covers only a small part of the phase space 
of interest to EAS. 

New p+C data at 12GeV/c and 158GeV/c, taken 
by the CERN experiments HARP g) at the PS ac- 



celerator and NA49 [7] at the SPS accelerator, are 
available now. At both beam momenta the sec- 
ondary particles (ir + ,ir~) are measured in a broad 
momentum region and up to large angles. At 
12GeV/c, the HARP data cover secondary momenta 
from 0.5 GeV/c to 8 GeV/c and an angular range from 
30mrad to 210mrad. The NA49 data are taken in the 
momenta range 0.85 GeV/c < p < 82.6 GeV/c and at 
angles up to 440 mrad 0. In addition there exist also 
several p+Be data sets (see references in (4)). How- 
ever, as simulations with hadronic interaction mod- 
els show, the particle production in p+air collisions is 
much more similar to p+C collisions than to p+Be re- 
actions because of the smaller difference in the atomic 
mass. 

4. Data from fixed target measurements 

4.1. NA49 p+C data at 158 GeV/c 

Pion spectra produced in p+C collisions with a 
beam momentum of 158GeV/c measured by NA49 
became available this summer (8). In Fig [3] rapidity 
spectra of secondary pions are compared with simu- 
lations done with QGSJET-01, QGSJET-II l9TT0ll and 
SIBYLL2. 1 HUH. SIBYLL2.1 and QGSJET-II pre- 
dictions show a reasonable agreement with the data, 
but QGSJET-01 overestimates the measured spectra 
by a factor ~ 1 .5 in the central rapidity region. Com- 
paring QGSJET-01 and QGSJET-II concerning the 
shape of the spectra, QGSJET-II predicts harder pion 
spectra than the older version. 

4.2. HARP p+C data at 12 GeV/c 

In these proceedings, we present for the first time 
momentum spectra of secondary pions in p+C colli- 
sions with a beam momentum of 12GeV/c measured 
by the HARP spectrometer (6). The thickness of the 
carbon target is equivalent to 5% nuclear interaction 
length (3.56 g/cm 2 ). For the selection of the sec- 
ondary particles (tt + , ir~) in the forward hemisphere 
four walls of drift chambers are used. The separa- 
tion of particle types is done with different detec- 
tor components: The time-of-flight measurement al- 
lows pion-kaon and pion-proton separation to be per- 
formed up to 3 GeV/c and beyond 5 GeV/c, respec- 
tively. Additional a Cherenkov detector is used to sep- 
arate pions from protons and kaons above 2.5 GeV/c 

na. 
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Figure 2. Coverage of the phase space regions of 
relevance to EAS (box histograms) by existing fixed 
target data using a proton beam and a carbon target 
(shaded/colored regions). Upper panel: total momen- 
tum of secondary pions vs. total momentum of proton 
projectiles. Lower panel: angle between beam and 
secondary particle momentum vs. beam momentum. 
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Figure 3. Rapidity spectra of pions produced in p+C 
collisions at 158GeV/c. The filled circles indicate 
the NA49 measurements [8|. The three lines show 
the results o f QGSJ ET-01 0, QGSJET-II l9ll0l and 
SIBYLL2.1 milia simulations. 

to 14% depending on the secondary momentum. Only 
for the highest momenta the statistical error of the ir~ 
amounts more than 20%. Generally the statistical er- 
rors of the 7r~ are some what larger than that for the 
7r + . The dominant contributions to the systematic er- 
ror are tertiary subtraction and momentum scale. 



In FiglH the preliminary tt + and tt~ momentum 
spectra are presented in six different angular bins 
starting from 0.03 rad up to 0.21 rad. The cross sec- 
tion of the 7T + production is consistently higher for all 
angular bins than the cross section of tt~ , especially 
in the higher momentum regions. This behaviour can 
be understood in terms of the leading particle effect 
which influences mainly positive pions. The error 
bars indicate the statistical and systematic errors. The 
size of both errors are of the same order of around 4 



5. Conclusions and outlook 

New p+C data from NA49 and HARP are avail- 
able now. They cover a broad range of phase space 
which is important for the EAS. However, only 20% 
of the grandmother particles in EAS are protons and 
most are pions. Currently an analysis of HARP ir+C 
data is in progress to help tuning hadronic interac- 
tion models for primary pions. Future data of the re- 
cently proposed upgraded MIPP experiment at Fermi- 
lab B14I15II of p+C, 7T+C and K+C collisions at 20, 
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Figure 4. Preliminary momentum spectra of secondary ir~ (filled circles) and tt + (open circles) mesons in 
p+C reactions at 12GeV/c measured with the HARP spectrometer at the PS accelerator at CERN. The six 
different panels show the spectra in different angular bins from 0.03 rad to 0.21 rad. 



60 and 120GeV/c would give important input for 
hadronic interaction models at higher energies. As 
a future project, an upgrade of the NA49 detector is 
planned ifTBI . This will allow the measurement of re- 
actions on p+C and ir+C at the energies of 30, 40, 50 
andl58GeV. 
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